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Abstract 


A standardized global data set of soil horizon thicknesses and 
textures (particle size distributions) has been compiled from the 

FAO/UNESCO Soil Map of the World, Vols . 2-10 (1971-81). This data set will 
be used by the improved ground hydrology parameterization (Abramopoulos et 
al., 1988) designed for the GISS GCM (Goddard Institute for Space Studies 
General Circulation Model) Model III . The data set specifies the top and 
bottom depths and the percent abundance of sand, silt, and clay of 
individual soil horizons in each of the 106 soil types cataloged for nine 
continental divisions. When combined with the World Soil Data File 
(Zobler, 1986), the result is a global data set of variations in physical 
properties throughout the soil profile. These properties are important in 
the determination of water storage in individual soil horizons and exchange 
of water with the lower atmosphere. The incorporation of this data set 
into the GISS GCM should improve model performance by including more 
realistic variability in land-surface properties. 
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Introduction 


As land- surface parameterizations in GCMs become more sophisticated, 

more detailed types of soil data are needed. Realistic models of water 
movement in the soil layer (e.g. Abramopoulos et al . , 1988) require 
information on variability of physical properties, including differences in 
the texture of soil horizons and thickness of these different horizons. In 
the global application of such land- surface parameterizations, the addition 
of geographic variability in the physical properties of soils, at the very 
least among the different continents , provides a more accurate description 
of surface characteristics that influence water movement in the soil 
profile. 

In the modified bucket parameterizations used in the GISS GCM Model 
II, soil water capacities are calculated as a function of vegetation type 
(Hansen et al . , 1983). More recently, global data sets of available soil 
water have been compiled using textural information from the upper 30 cm of 
the soil profile to assign water-holding attributes to the entire soil 
profile (Bouwman et al . , 1991) or the root zone (Henderson-Sellars et al., 
1986). In these studies, surficial physical properties are assumed 
representative of the entire soil profile and then estimates of soil water- 
holding capacity are assigned as a function of soil type and surficial 
texture for use in bucket-model calculations. In most soils, the textural 
class of the underlying horizons are significantly different from the 
surface texture due to pedogenic processes including degradation, 
aggradation, or translocation of materials within the soil profile (Buol et 
al., 1973). The bucket-model approach, however, is not appropriate when 
calculating water storage and movement throughout the complete soil profile 
as a function of hydraulic conductivity and matric potential. Furthermore, 
in a bucket model, infiltration can only be crudely estimated as a function 
of empirically determined parameters. These parameters, in turn, may 
unrealistically partition water at the atmosphere/soil interface into water 
lost as surface runoff versus water stored in the soil. 

We have generated a standardized global data set of physical soil 
properties consisting of soil horizon textures and thicknesses for 106 soil 
types based on information published by FAO/UNESCO (1971-81). The soil 
horizons are differentiated over the nine major continental divisions. 

This data set can be combined with Zobler's (1986; Figure 1) World Soil 
Data File to generate a global data set of physical properties that maybe 
used in simulating water movement in GCMs. 

The following report explains the decision-making process that went 
into compiling a standardized data set of texture and associated depth 
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information, the organization of the data set, and some of its strengths 

and weaknesses. We then compare three estimates for the potential storage 
of water in soil calculated using our data set with the field capacities 
used in the GISS GCM Model II. The three estimates are: 1) potential 

storage of water in the soil profile, 2) potential storage of water in the 
root zone, and 3) potential storage of water derived from soil texture. 

As investigators compiling a data set of physical soil properties to 
be used as prescribed surface characteristics in a GCM, we recognize that 
realistic values of available soil water can only be estimated as a 
combined function of climate, soil, and vegetation. Such research, though 
on-going, is beyond the scope of this report. 


Data and Methods 

The depth and textural (relative percent sand [2-0.05 mm], silt 
[0.05-0.002 mm], clay [<0.002 mm]) data for 106 soil types were taken from 
the Morphological, Chemical and Physical Properties Appendices in each of 
the nine volumes of the FAO-UNESCO Soil Map of the World (1971-81). These 
appendices listed data for selected soil profiles considered typical for 
each continental division. For some of the continental divisions the 
textural and depth data were not available for all of the soil types mapped 
by Zobler (1986). In other cases, the data were inconsistent or missing 
for one or more horizons of a soil type. 

We elected to include information from the A- , E- , B-, and C-horizons 
in our data set in recognition that these horizons contribute to the 
storage and movement of water within the soil profile. We did not include 
data from the litter (0-) horizon, from bedrock, or from parent material. 

We used data from the more complete description or the description with 
greater geographic coverage for those soil types with multiple descriptions 
of physical properties within a single continent. 


Decision Rules 

A set of decision rules was adopted to standardize the data set, to 
check the data for errors and to correct them, and to fill in missing data. 
Table 1 summarizes the adjustments that were made to data. 

1) Values were not included from horizons where the sum of percent sand, 
silt, and clay values was less than 50 or greater than 150 (presumably the 
result of either measurement or tabulation errors). 

2) A default basal depth of 3.6 meters was used when no bottom depth was 
specified for a soil profile. The default depth of 3.6 meters was selected 
to allow realistic simulation of dynamic hydrology. 
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3) An average depth was calculated in cases when a depth range was given 
or the top and bottom depths of contiguous horizons were not the same. 
Depths reported in inches were converted to metric. 

4) To use the particle size information for calculation of hydraulic 

conductivity and matric potential, the percent sand, silt, and clay data 
had to sum to 100 percent. To normalize the data to sum to 100 percent, we 
adjusted the percent clay values, rather than renormalizing all three size 
classes, because of the greater potential error in the measurement of the 
clay size class. As long as the summed percent ranged between 80 and 120 
percent, the appropriate amount was subtracted/added from the clay 
fraction. In one case, when the summed percent was more than 120 percent 
or less than 80 percent, however, the qualitative descriptions of the soil 
horizons were used as guidelines to determine the proportions of each size 
class . 

5) Interpolation of values as the average of data from bracketing horizons 
was used to fill in missing data for individual horizons. We elected to be 
conservative and presuppose continuity of trends within the soil profile 
rather than prescribe data on the basis of qualitative descriptions. We 
recognize that assigning values to a B-horizon as the intermediate 
properties of A- and C-horizons is not representative of many soils but we 
believe that this is more favorable to defining values arbitrarily. 

5a) When qualitative descriptions suggested that adjacent horizons were 
very similar (e.g., similar horizon nomenclature: Bgl, Bg2, Bg3, ...), data 
from the adjacent horizon were substituted for the missing data. 

6) Extrapolation of values assuming a linear rate of change was used to 
replace missing values in either the uppermost or lowermost horizon. Once 
again this is a conservative choice because we assumed continuity of trends 
between adjacent soil horizons. 

7) Adjustments/Corrections were made for obvious tabulation errors and 
large discrepancies between the qualitative descriptions and percent sand, 
silt, and clay values. For many of the North Central Asia soils and a few 
soils from the other continental divisions, the particle size analyses were 
reported in nonstandard international particle size ranges. These data 
were converted to standard size ranges (sand > 0.05 mm; silt 0.05-0.002 mm; 
clay < 0.002 mm) by assuming uniform distributions within each reported 
size range and partitioning the data proportionally . 


Hissing Data 

A number of soil types from each continent were completely missing 
data. To fill in the textural or depth data for these cases, we 
substituted data using the same soil type from a different continental 
division. The substitution order for each continent was made recognizing 
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similarities in geology and modern climates. Table 2a lists the hierarchy 
used to substitute data from a different continent. When data for a soil 
type were absent from all the continents, data were substituted from an 
adjacent soil type with similar descriptive characteristics. Table 2b 

lists these soil types and their replacement soil types. Lithosols were 
assigned a total thickness of 10 cm and the particle size distributions of 
the C-horizon of a Eutric Regosol. Rankers were assigned a total thickness 
of 25 cm and the particle size distributions of the A-horizon of a Humic 
Cambisol. Histosols by definition are organic material and cannot be 
described in terms of relative percentage of sand, silt, and clay particle 
sizes. We assigned Histosols an average thickness of 360 cm and arbitrary 
distribution of particle sizes summing to 100 percent but not used in any 
subsequent calculations . 

The frequency of occurrence of each of the 106 soil types for each of 
the nine continental divisions and missing data substitutions are listed in 
Table 3. In most cases, the potential error associated with substituting 
data for soil types with missing information is not large. The frequency 
of occurrence for these soils in Zobler's World Soil Data File rarely 
exceeds 50 lxl grid cells, excluding Lithosols, Rankers, and Histosols. 

The most abundant soil type substituted for was Haplic Yermosol, with 615 
lxl grid cells. Among the continental divisions, the absence of soil 
properties data was most widespread for Mexico/Central America followed by 
North-Central Asia. We did not substitute data for the missing values of 
two minor soil types, Ferric Podzol and Gelic Planosol, since these were 
not present in the Zobler World Soil Data File. 

Soil Profile Thickness 

Maximum soil depths in the data set are shown in Table 5 and Figure 
2. In many cases, the soil profile thicknesses represent minimum possible 
values because profile descriptions do not always extend to subsurface 
bedrock. The soil thicknesses range from 10 cm for Lithosol to 800 cm for 
Distric Nitosol in Africa. The spatial distribution of soil profile 
thickness can be summarized as thickest in the well-developed soils of 
tropical low latitudes and thinnest in the poorly developed soils of high 
latitudes. The soil profiles are thin in mountainous regions such as the 
Himalayas and Andes and are thick in mid- latitude peatlands such as those 
found in northern Europe and North America. 

Potential Water Storage Calculations 

Using the data set, we have calculated three different estimates of 
potential storage of water in soil. These estimates can be used as proxies 
for the amount of water in the soil layer available to plants for 
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evaporation. They were all estimated as the sum of the amount of water in 
a soil column with a 1x1 cm cross section. 

1) An estimate of potential storage of water in the soil profile was 
calculated. The relative saturation capacities for the different particle - 

size classes and for peat (Histosols) listed in Table 4a. 

2) An estimate of potential storage of water in the root zone was also 
calculated. The root-zone thickness was derived using a) information from 
the simplified 8-type version of the Matthews (1983) 22-type global 
vegetation data set (Matthews, 1984; Hansen et al . , 1983) and b) estimates 
of maximum root-zone thicknesses for these eight types (Table 4b; 
Rosenzweig, unpub.). When the soil profile thickness is larger than the 
maximum root-zone thickness, root-zone thickness is limiting. 

3) The potential storage of water derived from soil texture was calculated 
as a function of the textural class (Table 4c) of each horizon within the 
soil profile. The textural class of each horizon was estimated from the 
relative amounts of sand, silt, and clay. 

These estimates of potential storage of water in soil are compared 
with the global distribution of the vegetation- dependent water-holding 
capacity prescribed for the GISS GCM Model II (see Table 6 in Hansen et 
al., 1983). The water-holding capacities in Model II range from 20 mm for 
desert vegetation to 650 for cm rainforest vegetation. The spatial 
patterns of this estimate closely correspond to the mapped patterns in 
Matthews (1983) 22 -type global vegetation data set. Low values are 
concentrated in subtropical and polar desert regions, whereas high values 
are associated with low latitude rainforests. 


Derived Properties 
Potential Storage of Hater in the Soil Profile 

The potential storage of water at 100 percent saturation represents 
the maximum amount of water that the soil profile can possibly hold (Table 
6; Figure 3a). The values for potential storage of water for the entire 
soil profile range from 42 mm for Lithosol to 4432 mm for Distric Nitosol. 
The geographic distribution of this estimate shows that areas of high 
values closely correspond to areas of thick soil profiles with high clay 
content resulting from greater soil development; likewise, the areas of low 
values are associated with thin soil profiles with low clay content. 

Distinct features include low values for high latitude Lithosols in 
northeastern Asia, northern North America, mountainous regions, and central 
desert areas. Large values are located along the equator in South America, 
Africa, and south Asia, and in midlatitude peatlands of northern Europe and 
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North America. The major desert regions are not resolved with this 
estimate because of the relatively large soil horizon thicknesses, 
underscoring that potential storage of water in the soil profile is the 
maximum amount of water that can be stored throughout the soil profile. 

Overall, the potential storage of water in the entire soil profile over 
estimates the amount of soil water available for evapotranspiration. Soils 
are rarely completely saturated throughout the profile. 


Potential Storage of Water in the Root Zone 

The potential storage of water in the root zone is shown in Figure 
3b. This measure indirectly includes climate information as a function of 
the vegetation data because vegetation coverage and associated root- zone 
thickness reflect climatic moisture gradients. For example, maximum 
rooting depths in desert regions are determined by water supply. The 
values of the potential storage of water in the root zone range from less 
than 2 mm for desert type vegetation associated with Lithosols, Arenosols, 
and Xerosols to as large as 1700 mm for Woodland, Evergreen, and Deciduous 
vegetation associated with Histosols. 

For the most part, the mapped pattern of this measure resembles the 
Matthews (1983) 22-type global vegetation data set. Low values associated 
with major desert areas are well defined in north Africa and central Asia. 
The low values farther north correspond to high latitude Lithosols . High 
values in the moderate -to -high latitudes are often associated with 
forested, deep organic-rich soils. The low values in tropical rainforest 
areas of northeastern South America and west- central Africa derive from the 
moderately shallow maximum rooting depth used for rain forest vegetation. 


Potential Storage of Water Derived from Soil Texture 

A commonly used estimate of the amount of soil moisture available for 
evapotranspiration is calculated as a function of texture (Table 7; Figure 
3c) . This estimate is defined as the amount of water released between in 
situ field capacity and the soil wilting point (usually measured as the 
difference in water content at soil matric potentials of -0.03 MPa and 
-0.15 MPa; Soil Science Society of America, 1987). The values for 
potential storage of water derived from soil texture range from 17 mm for 
Lithosol to 2160 mm for Histosols. In general, the geographic distribution 
of this estimate shows that areas of high values also correspond to areas 
of organic soils or thick soil profiles with high clay content resulting 
from greater soil development. 


6 


Comparisons with the GISS GCM Model II Water-Holding Capacities 

Spatial variation in the Model II water-holding capacities (Figure 
3d), though consistently lower, resembles the estimate of potential storage 
of water in the root zone (Figure 3b) . The similarities in geographic 

variability result from both the Model II and the root- zone estimates being 
derived from Matthews (1983) vegetation data set. The maximum possible 
water -holding capacity used in Model II of the GISS GCM (650 cm) is 
significantly smaller than the maximum possible values for either of our 
three estimates of potential storage of water in soil (1700-4432 cm) . The 
mapped values of water-holding capacities used in Model II are concentrated 
at the low end of the scale and show much broader and smoother patterns of 
spatial variability than the three estimates we calculated. 


Discussion 

The task of generalizing soil properties for GCMs requires 
transforming very heterogeneous high resolution data into representative 
homogeneous information at a much lower resolution. Variability in soil 
physical properties (e.g. texture and thickness) is scale - dependent . While 
there remains a great deal of variability at the subcontinental scale not 
captured in our data set, we believe characterizing these properties at the 
continental scale is a necessary first step to improve current 
specifications of soil in GCMs. 

A number of less -than- satisfactory GCM results have been attributed 
to inadequate ground hydrology parameterizations and the values used for 
soil water-holding capacity (Rind, 1988; Rind et al . , 1990; Delworth and 
Manabe, 1988; Kellogg and Zhao, 1988). For example, Rind (1988) has 
suggested that primitive ground hydrology parameterizations in the GISS GCM 
have resulted in "too much soil moisture and rainfall" for Canada and the 
northern USA. The data set we have compiled will be used in an updated 
version of the improved ground hydrology parameterization that is being 
incorporated into the current version of the NASA/GISS GCM (Abramopoulos et 
al., 1988). We expect that GISS GCM experiments using the new ground 
hydrology scheme and the more realistic soil characteristics from our data 
set will provide insights into climate model sensitivity to energy/moisture 
fluxes between the land and atmosphere and lead to improved model 
performance . 

Examination of the differences among our three potential storage of 
water in soil estimates demonstrates the problems that one may encounter 
portraying soil moisture as a static, empirically-derived estimate. Soil 
moisture is an extremely interactive measure reflecting the collective 
influence of climate, soil properties, and vegetation. Focusing on our 
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results for northern Africa helps to illustrate this point. When the 
entire soil profile is considered as a possible reservoir, the soil profile 
has a relatively large storage capacity. In contrast, when only the 
maximum rooting depth is considered as a possible reservoir, the soil 
profile has an extremely small storage capacity. The soil texture-based 
water- storage estimate represents an intermediate alternative with 
moderately low storage capacity. The GISS Model II water-holding 
capacities for the region, also derived from vegetation data, are similar 
to the values based upon rooting depth. Under current dry conditions, the 
very small capacity seems intuitively more correct; however, if climate 
change resulted in regionally wetter conditions, the intermediate or the 
large capacity would be more appropriate. 

In addition to describing our data set of physical property of soils, 
we have presented estimates for soil moisture as the potential storage of 
water in the soil profile, the potential storage of water in the root zone, 
and the potential storage of water derived from soil texture. 

Nevertheless, we feel it is important to caution that realistic estimates 
of soil moisture can only be calculated using an interactive climate 
forcing of moisture supply and demand. 
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Appendix 1: Description of the Computer Files 


The particle size data set has been archived at the National Space 
Science Data Center, NASA/Goddard Space Flight Center, Greenbelt, MD 

20771. The data has been stored in free format as four 106x10x15 
dimensioned real*4 arrays: depth, sand, silt, and clay. The first 
dimension (106) corresponds to the sequence number of the soil types in 
Zobler's (1986) World Soil Data File. The second dimension (10) 
corresponds to the volume numbers of the nine major continental divisions 
in the FAO/UNESCO Soil Map of the World, Vols. 2-10 (1971-81). The third 
dimension (15) corresponds to the individual horizons with data for each 
soil type from the Morphological, Chemical and Physical Properties Appendix 
in each of the nine volumes of the FAO/UNESCO Soil Map of the World (1971- 
81). The data in the sand, silt, and clay arrays are stored as 
proportional values for each soil horizon. The arbitrary particle size 
distribution summing to 100 percent included for Histosols (entries 61-63 
in the first dimension of each array) should not be used. Instead, values 
reflecting the physical properties of organic soils and appropriate for 
specific research objectives should be inserted. 

The data in the depth array are scaled in meters with the first value 
being 0 m depth for each soil type and the subsequent values the contact 
depths of contiguous horizons. By definition, the depth array contains one 
extra value for the third dimension corresponding to the bottom depth of 
the lowest horizon for each soil type. Within the data set, no soil type 
had more than 14 soil horizons. In cases when the number of horizons in a 
soil type was less than 14, we used -1.0 values to flag the end of the 
record for each soil type. For example, a soil type with 10 horizons has 
10 data entries in the sand, silt, and clay arrays, 11 data entries for the 
depth array, and -1.0 values for entries 11 - 15 in each array (entries 12 
- 15 for the depth array) . 
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Figure Captions 


Figure 1. Global distribution of 26 major soil units mapped by Zobler 
( 1986 ) . 

Figure 2. Global distribution of soil profile thickness based on maximum 
soil depth of each soil type. 

Figure 3. Global distribution of a- the potential storage of water in the 
soil profile, b- the potential storage of water in the root zone, c- 
the potential storage of water derived from soil texture, and d- 
water-holding capacity prescribed for the GISS GCM Model II. 
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depth data for one or more horizons in profile corrected/adjusted 

clay particle size data for one or more horizons in profile normalized to sum to 100 percent 

particle size data for one or more horizons in profile interpolated 

particle size data for one or more horizons in profile extrapolated 

particle size data for one or more horizon in profile corrected/adjusted 


Table 2a. Hierarchy for substituting for missing data using information from a different continent 
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Table 3. Summary of frequency of occurence of each soil type in Zobler World Soil File with numeric and 
letter codes in parentheses flagging the type of missing data substitution. 
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10 Australia/South Asia 

Two letter codes refer to Column 2 within table 


Sand 0.394 (Zobler, unpub. 
Silt 0.537 (Zobler, unpub. 
Clay 0.577 (Zobler, unpub. 
Peat 0.85 (Webb, 1990) 


w 

<D 

W 

W 

03 


Q) 

pH 

X 

03 

H 



X 

£ 

a> 

X 

a) 

M-t 


Sh 

O 


a> 

x 

o3 

£ 


a> 


% 

o 

X 

§ 

0 

1 


> 

•H 

4-> 

Cti 


a 


xi 



20 






Table 5. Soil thickness for 106 soil types, 9 continents (cm) 
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Table 6. Potential storage of water in the soil profile for the 106 soil types, 9 continents (mm) 
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Figure 1 . Soil Profile Thickness. 
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Figure 2. Major Soil Units 
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A) POTENTIAL STORAGE OF WATER IN SOIL PROFILE 


B) POTENTIAL STORAGE OF 


WATER IN ROOT 


ZONE 


Figure 3a and 3b. 


32 






FOLDOUT F RA 1 ^ £QLDOUT FRAME ^ 

□OUT fWM 


r 




foldchjt frame foldout framb 




Photoheliograph :'hase C and P Program 
Work Breakdown Structure 



Figure 3c and 3d. 
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